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By using the dual vortex method (DVM), we study some phases and phase transitions in an 
extended boson Hubbard model slightly away from 1/3 ( 2/3 ) filling on frustrated lattices such 
as triangular and kagome lattice. We develop systematically a simple and effective way to use the 
vortex degree of freedoms on dual lattices to characterize the symmetry breaking patterns of the 
boson insulating states in the direct lattices. In the triangular lattice at 1/3, we find a X-CDW, 
a stripe CDW phase which was found perviously by a density operator formalism (DOF). Most 
importantly, we also find a new CDW-VB phase which has both local CDW and local VB orders, in 
sharp contrast to a bubble CDW phase found previously by the DOF. In the Kagome lattice at 1/3, 
we find a VBS phase and a 6 fold-CDW phase. Most importantly, we also find the CDW-VB phase 
which has both local CDW and local VB orders which was found in previous QMC simulations. We 
also study several other phases which are not found by the DVM. By analyzing carefully the saddle 
point structures of the dual gauge fields in the translational symmetry breaking side, we find that the 
CDW to SF transition in the triangular is likely to be a strong first order one, while the CDW+VB 
to the SF transition in Kagome lattice is likely to be a weakly first order one. This findings by 
the DVM can give some explanations to some puzzles observed in the QMC in both lattices. Then 
by pushing the effective actions slightly away from the commensurate filling / = 1/3, we classified 
all the possible types of supersolids and analyze their stability conditions. In a triangular lattice, 
there are X-CDW supersolid, stripe CDW supersolid, but absence of any valence bond supersolid ( 
VB-SS ). There are also a new kind of supersolid: CDW-VB supersolid. In a Kagome lattice, there 
are 6 fold-CDW supersolid, stripe CDW supersolid, but absence of any valence bond supersolid ( 
VB-SS ). There are also a new kind of supersolid: CDW-VB supersolid. We show that independent 
of the types of the SS, the quantum phase transitions from solids to supersolids driven by a chemical 
potential are in the same universality class as that from a Mott insulator to a superfluid, therefore 
have exact exponents z = 2, v — 1/2, r\ = ( with logarithmic corrections ). Implications on QMC 
simulations with both nearest neighbor and next nearest neighbor interactions in both lattices are 
given. Some possible problems of the DOF in identifying the insulating phases are also pointed out. 



I. INTRODUCTION 

The EBHM with various kinds of interactions, at 
various kinds lattices ( bi-partisan or frustrated ) 
at various kinds of filling factors ( commensurate 
or in-commensurate ) is described by the following 
Hamiltoniar>ir— : 

H = -t (b\bj + h.c.) - ^y^i + - 1) 

<ij> i i 

+ Vi UiUj + V 2 ^2 riin k H (1) 

<ij> <<ik» 

where rii — b\bi is the boson density, t is the hopping 
amplitude, U,V\,V% are onsite, nearest neighbor (nn) 
and next nearest neighbor (nnn) interactions between 
the bosons. The • ■ ■ may include further neighbor in- 
teractions and possible ring-exchange interactions. A su- 
persolid in EqnJT]is defined as to have both off-diagonal 
long range order < b{ >^ and diagonal charge den- 
sity wave in the boson density n, which break the lattice 
symmetry. In frustrated lattices, the sign of the hopping 
t in EqnQ] makes the crucial difference- This difference 
leads to many new and interesting physics in frustrated 



lattices which are quite different than those in bipartite 
lattices studied 

In cold atom experiments^, the on-site interaction U 
can be tuned by the Feshbach resonance. Various kinds 
of optical lattices can be realized by suitably choosing 
the geometry of the laser beams forming the optical lat- 
tices. For example, using three coplanar beams of equal 
intensity having the three vectors making a 120° angle 
with each other—, the potential wells have their minima 
in a honeycomb or a triangular lattice. Four beams travel 
along the three fold symmetry axes of a regular tetrahe- 
dron, the potential wells have their minima in a body- 
centered-cubic lattice^. The authors in£ proposed to cre- 
ate a kagome optical lattice using superlattice techniques. 
There are many possible ways to generate longer range 
interaction Vi , Vi , ■ ... of ultra-cold atoms loaded in optical 
lattices (1) Very exciting perspectives have been opened 
by recent experiments on cooling and trapping of 52 CY 
atomsiS and fermionic polar molecules A0 K + 87 Rb^. The 
bosonic polar molecules 39 K + 87 Kb is also expected in 
near future. Being magnetically or electrically polar- 
ized, the 52 Cr atoms or polar molecules interact with 
each other via long-rang anisotropic dipole-dipole inter- 
actions. Loading the 52 Cr or the polar molecules on a 2d 
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optical lattice with the dipole moments perpendicular to 
the trapping plane can be mapped to Eqnfl] with long- 
range repulsive interactions ~ p 2 /r 3 where p is the dipole 
moment. The CDW supersolid phases described by the 
dual vortex method (DVM) ir£ was numerically found 
to be stable in large parameter regimes in this system^. 
Possible techniques to generate long-range interactions 
in a gas of groundstate alkali atoms by weakly admixing 
excited Rydberg states^ with laser light was proposed 
inM. Very recently, the cold bosons or fermions in opti- 
cal lattice can interact with a very long range interaction 
due to the exchange of cavity photons 1 —. The generation 
of the ring exchange interaction has been discussed iitiS. 

From a spin wave 1 / S expansion, the authors ir>2 found 
that a SS state is more robust in a triangular lattice 
with only t and V\ terms in EqnQ] and a SS state with 
\[2> x v3 pattern is stable even at half filling / = 1/2. 
Indeed, this discovery was confirmed by several recent 
QMC simulations in the EBHM of hard core bosons with 
the nearest neighbor (NN) interaction^—. However, 
due to too strong quantum fluctuations, the authors ir*2 
found that the spin wave expansion does not work any- 
more in Kagome lattice. So a different analytical method 
is needed to study possible supersolid phases in a Kagome 
lattice. 

Recently, the EBHM in a square lattice at generic com- 
mensurate filling factors / = p/q ( p, q are relative prime 
numbers ) were systematically studied in 1 - by a dual vor- 
tex method ( DVM ). The DVM method was applied by 
one of the authors to study the EBHM on bipartite lat- 
tices such as honeycome and square lattice at half and 
near half fillings and its possible experimental applica- 
tions in cold atoms, adatom absorptions on substrates 
and possible Cooper-pair supersolids in high tempera- 
ture superconductors. It was pointed out ir£ that the 
DVM developed in^ holds in the superfluid ( SF ) and 
the valence bond solid ( VBS ) side where the saddle 
point of the dual gauge field can be taken as uniform, 
but cares are needed to apply the DVM developed ini in 
the CDW side where the saddle point of the dual gauge 
field can not be taken as uniform anymore. When study- 
ing phases and phase traditions slightly away from 1/2 
filling, special care is needed to choose a correct sad- 
dle point of the dual gauge field in the CDW side, so 
a different effective action is needed in the CDW side 
to make the theory self-consistent^. In£, (1) by extend- 
ing the DVM explicitly to the lattice symmetry breaking 
side by choosing the corresponding self-consistent saddle 
points of the dual gauge field (2) by pushing the DVM 
to slightly away from commensurate filling factors, the 
author mapped out the global phase diagram at T = 
of the chemical potential versus the ratio of the kinetic 
energy over the interaction. It was found that in the 
insulating side, different kinds of supersolids are generic 
possible states slightly away from half filling. A novel 
kind of supersolid called valence bond supersolid ( VB- 
SS ) was proposed. It was also shown that the transition 
from a CDW ( Valence bond ) solid to a CDW ( valence 



bond ) supersolid driven by the chemical potential is in 
the same universality class as the Mott to the SF with 
the exact exponents z = 2,v = 1/2,7] = subject to 
a logarithmic correction. It was also suggested that the 
density order parameter constructed ir*i to characterize 
the symmetry breaking patterns in the insulating states 
are also not physically transparent. It is very hard to 
generalize the density operator to a honeycomb lattice 
and other lattices. The subsequent QMC simulations on 
the EBHM of hard and soft core bosons with the near- 
est neighbor (NN) interaction on a honeycomb lattice^ 
indeed found a stable X-supersolid phase in the soft core 
case slightly above the half filling. 

The DVM method was extended to study the EBHM in 
a triangular lattice at 1/3 and 1/2 fillings ir*2i. Starting 
from the SF side, the authors derived the effective action 
Eqn|3] to describe quantum phase transitions from the 
SF to some insulating states with several kinds of lattice 
symmetry breaking patterns. By generalizing the density 
operators in the square lattice constructed ini to those 
in a triangular lattice, they identified an X-solid phase at 
/ = 1/3 in the Ising limit of the effective action, a stripe 
solid phase and a bubble solid phase in two possible easy 
plane limits. They suggested that the transition from 
the SF to the stripe solid could be second order through 
the so-called de-confined quantum critical points. Then 
subsequent QMC simulations on the EBHM Eqn|T2] of 
hard core bosons with a next nearest neighbor (NNN) 
interaction ir>22 indeed found the stripe solid phase as a 
ground state at / = 1/3 in some parameter regimes. The 
QMC in32i seems also to find a meta-stable bubble phase 
which has higher energy than the stripe phase in the 
same parameter regimes. However, they found the tran- 
sition from the SF to the stripe solid phase at / = 1/3 is 
a strong first order transition. Later, the DVM was ap- 
plied to study the EBHM on a Kagome lattice in2si. From 
the Magnetic space group (MSG), starting from the SF 
side, the authors derived the effective action EqnJ3] to 
describe quantum phase transitions from the SF to some 
insulating states with some lattice symmetry breaking 
patterns. Realizing it is not convenient to generalize the 
density operators in the square lattice constructed ini 
to the Kagome lattice, they tied to use the vortex fields 
on the dual dice lattice to identify possible symmetry 
breaking patterns in the insulating side. Unfortunately, 
because the vortex fields are gauge dependent, so may 
not be used to characterize symmetry breaking patterns 
effectively. There are also two QMC simulations on the 
EBHM of hard core bosons with nearest neighbor (NN) 
interaction in a Kagome lattice2&22. Both groups found 
that the solid state at / = 1/3 has both CDW and VB or- 
der. The VB order corresponds to boson hopping around 
a hexagon ( see Fig. 13b). However, they did not see any 
stable supersolid phase away from the 1/3 filling in the 
hard core case. Although the QMC in Refj^l suggested 
that the the SF to the CDW+VB solid transition is a 2nd 
order transition through the deconfined quantum critical 
point^, the QMC in Ref*2!l concluded a weakly 1st order 
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transition from a double-peaked histogram of the boson 
kinetic energy. 

In this paper, we will use the DVM developed ini^ to 
study quantum phases, especially various kinds of insu- 
lating phases and phase transitions in the two most com- 
mon frustrated lattices such as triangular and kagome 
lattices at and slightly away 1/3 and 2/3 fillings. We con- 
trast different symmetry breaking patterns of the insulat- 
ing phases in both Ising and Easy plane limits in the two 
lattices. We will also compare our results with some pre- 
vious Quantum Monte-Carlo (QMC) results in the two 
lattice a 23 i 26 and also give important implications to pos- 
sible future QMC simulations on the EBHM EqnQ]with 
both hard and soft cores, with both the NN V\ and the 
NNN V 2 interactions. We develop a systematic way to de- 
termine the symmetry breaking patterns of the insulating 
states in terms of the vortex degree of freedoms only at 
the dual lattice points. These vortex degree of freedoms 
are the gauge invariant physical vortex densities, the ki- 
netic energies and vortex currents defined in Eqn l7l8l In 
the triangular lattice, we found the X-solid phase in the 
Ising limit v > 0, the stripe phase in one of the easy 
plane limit v < 0, w < 0, both of which were found be- 
fore by the density operator formalism in^L. However, we 
identify a novel phase with both CDW and VB orders 
shown in Fig. 7 in another easy plane limit v < 0,w > 0. 
This finding differs from the bubble solid phase shown in 
Fig. 25 in the same easy plane limit found by the density 
operator formalism in 21 . This disagreement calls for the 
reexamination of the correctness and effectiveness of the 
density operator formalism in terms of the dual vortex 
degree of freedoms in general lattices at general filling 
factors. The quantum phases and phase transitions in 
a Kagome lattice in both Ising and easy plane limit are 
dramatically different from those in bipartite^^ and tri- 
angular lattice studied previously 2 ^- and in the Sec. II. So 
far, the density operator formalism in£i2i has not been 
generalized to the Kagome lattice However, by using our 
method, we identify a triangular valence bond (TVB) 
solid phase in the Ising limit v > 0, a 6-fold CDW phase 
in one of the easy plane limit v < 0, w < 0, the CDW- 
VB phase first discovered in the QMC i n 26 ' 27 in the other 
easy plane limit v < 0, w > 0. In fact, this 6-fold CDW 
phase has the same symmetry breaking patterns as the 
CDW-VB phase, however, the crucial difference is that 
the former has a local CDW order, while the latter has 
a local VBS order. In the Ising limit, the SF to the 
triangle valence bond order (TVB) is a weak first or- 
der one. For the first time, the Ising limit of the same 
action can be used to describe a VBS phase in a con- 
crete model. In the easy-plane limit, we find that the 
6 fold CDW to the SF transition, also the CDW+VB 
to the SF transition at commensurate filling / = 1/3 in 
Kagome lattice are a weakly first order one. This finding 
resolves some outstanding puzzles observed in the QMC 
in a kagome lattic o 26 ' 27 . The firm identification of this 
CDW-VB phase in the easy-plane limit v < 0, w > is 
important, because one is sure that the effective action 



Eqn|T3] in terms of the two vortex fields 4>o,(f>i in the 
easy-plane limit v < 0, w > indeed describes the SF to 
the CDW-VB transition studied by the QMC 2 ^ 2 - 7 - in a 
simple microscopic EBHM Eqn[TT] in a Kagome lattice. 
This EBHM EqnQT] is very simple and has no ring ex- 
change interaction which is usually needed to stabilize a 
VBS order in a bipartite lattice. The results ir>2& showed 
that this transition is a very weak first order one instead 
of a second order one through the DCQCP. We also study 
briefly the excitation spectra in these insulating phases. 
It is interesting to note that one can achieve different 
CDW+VB states in the easy plane limit v < 0, w > in 
both triangular and Kagome lattices from the DVM. The 
method developed here should be very general and can 
be used to characterize uniquely the symmetry breaking 
patterns in any lattices at any filling factors by using the 
vortex degree of freedoms only at the corresponding dual 
lattice points. 

Then we study possible quantum phases and quantum 
phase transitions slightly away from 1/3(2/3) fillings in 
the two frustrated lattices. When pushing the effective 
action slightly away from the commensurate fillings in- 
side various insulating phases, we find that one has to 
choose the corresponding saddle point structure of the 
dual gauge field carefully and self-consistently to con- 
struct the effective actions away from the commensurate 
fillings on insulating side. We classified all the possible 
types of supersolids (SS) and analyze their stability con- 
ditions. In a triangular lattice, there are X-CDW SS and 
stripe CDW SS which are stable in both hard core and 
soft core cases, but absence of any Valence bond super- 
solid ( VB-SS ). There are also a new kind of supersolid: 
CDW-VB supersolid which maybe stable only in the soft 
core case. This novel kind of supersolid has CDW, VBS 
and superfluid order. The supersolid in the triangular 
lattice near 1/2 filling can be considered as doping the 
adjacent solid at 1/3 filling by interstitials ( SS-i in Fig. 9) 
or as doping the adjacent solid at 2/3 filling by vacan- 
cies ( SS-v in Fig. 9). The supersolid at exactly 1/2 is 
the coexistence of SS-i and SS-v with any possible ratio. 
There are only two kinds of supersolids: vacancy type or 
interstitial type. There is no other kinds of supersolids. 
The 1/2 filling at triangular lattice may not be a " com- 
mensurate " filling as thought previously. Only 1/3(2/3) 
are commensurate fillings. Therefore the SS in triangular 
lattice has the same origin as that in the square lattice 
discussed in£. In a Kagome lattice in the Ising limit, 
there can only be a direct first-order transition from the 
triangular valence bond (TVB) to the SF, so there is no 
immediate TVB supersolid intervening between the TVB 
and the SF. However, in the Easy-plane limit, different 
kinds of supersolids are generic states slightly away from 
1/3 ( 2/3 ) filling. In addition to a 6-fold CDW super- 
solid, the stripe supersolid ( stripc-SS), we also find a new 
kind of supersolid: CDW-VB supersolid. The stripe-SS 
should be stable even in the hard core limit, while the 
stripe-SS and CDW- VB-SS maybe stable only in the soft 
core limit. The absence of the VB-SS in frustrated lat- 
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tices is in sharp contrast to that in bipartite lattices dis- 
covered ini£. The quantum phase transitions from solids 
to the adjacent supersolids ( if they exist ) driven by a 
chemical potential are in the same universality class as 
that from a Mott insulator to a superfluid, therefore have 
exact exponents z — 2, v — 1/2, r\ = ( with logarithmic 
corrections ). The superfluid density in the SS scales as 
p s ~ \p- l/3\( d+z - 2 > = |/5-l/3| = Sf with logarithmic 
corrections. In the anisotropic stripe SS case, the super- 
fluid densities along different directions should scale the 
same way with different coefficients. We also made im- 
plications on the available and future QMC simulations 
on EBHM with both NN and NNN interactions in both 
lattices. 

The rest of the paper are organized as follows. In the 
section II, we study the EBHM in a triangular lattice. 
We first outline the method to characterize the symmetry 
breaking patterns in the insulating side using the vortex 
degree of freedoms only at dual lattice sides. Then we 
identify the X-CDW phase in the Ising limit in Sec. II- 
A, the stripe CDW phase in Sect.II-B-1, the new and 
the most interesting CDW-VB solid phase in Sect.II-B- 
2, the VBS phase in Sect.II-C, then compare with QMC 
in Sec.II-D. We also construct effective actions to study 
possible supersolid phases and their stabilities slightly 
away from the 1/3 filling. In the section III, we study 
the EBHM in a Kagome lattice. We identify the VBS 
phase in the Ising limit in Sec.III-A, the most interesting 
CDW-VB phase in Sect.III-B-1, the 6 fold CDW phase 
in Sect.III-B-2, the stripe CDW phase in Sect.III-C, then 
compare with QMC in Sec.III-D. We also construct effec- 
tive actions to study possible supersolid phases and their 
stabilities slightly away from the 1/3 filling. We summa- 
rize our results and give some perspectives on compar- 
ing the DVM and the QMC simulations in the conclud- 
ing Sect. IV. In the appendix A and B, we use the same 
method to identify insulating phases in square lattice and 
honeycomb lattice respectively, also the excitation spec- 
tra in all these phases. In appendix C, we compare our 
method with the density wave operator formalism used 
and commented on possible weakness of the formal- 
ism. In appendix C, we list specific vortex fields values 
used in Sect. II-B-2 for the CDW-VB phase. A very 
short version in the part of Kagome lattice appeared in 
a unpublished preprint 28 . 



II. SOLID AND SUPERSOLID NEAR 1/3(2/3) 
ON TRIANGULAR LATTICE 

Triangular lattice is the simplest frustrated lattice ( 
Fig. la). From spin wave expansion, the authors in 2 found 
that a SS state is more robust in a triangular lattice 
with only t and V± terms in Eqn[T] and a SS state with 
\/3 x V3 pattern is stable even at half filling q = 2. 
Indeed, this discovery was confirmed by several recent 
QMC simulations^^. 

We will investigate the phases and quantum phase 



transitions when slightly away from the 1/3 filling. 
Because of the P-H symmetry at U — oo, the re- 
sults are equally valid near 2/3 filling. As shown in 
the Table 1 in 3 , for q = 3, there are 3 minima at 
(0,0), (0,27r/3), (0,47r/3). Let's label the three eigen- 
modes at the three minima as I — 1,2,3. The general 
effective action in terms of the three modes invariant un- 
der all the MSG transformations upto sixth order terms 
was written downSI. The effective action can be simpli- 
fied in the permutativc representation basis <pi,l — 1, 2, 3 
given by: 

1 -2* 

Co = -^{4>o+e^ 0i + 2 ) 

1 .2, 

6 = -y=(e l 3 00 + ^1 + 02) 

6 - -±=(0o + 0i+e^0 2 ) (2) 

In the superfluid side, moving slightly away from the 
1/3 filling f — 1/3 corresponds to adding a small mean 
dual magnetic field 5f = f — 1/3 in the action derived 
in 21 . Upto sixth order, the action is Csf = Cq + Ci+ C-2- 

£o = ^|(a M -a M )0 ; | 2 +r|0 / | 2 + (e^ A a,AA-27r ( 5/V) 2 /4 
I 

& = U (|0 O | 2 + |0l| 2 + |02| 2 ) 2 -«[(|0o| 2 -|0l| 2 ) 2 

+ (|0i| 2 -|0 2 | 2 ) 2 + (|0 2 | 2 -|0o| 2 ) 2 ] 
C 2 = wMoM 3 + (0i 2 ) 3 + (02>o) 3 + h.c] (3) 

where is a non-compact U(l) gauge field. 

Because the duality transformation is a non-local 
transformation, the relations between the phenomenolog- 
ical parameters in Eqn|3]and the microscopic parameters 
in Eqn[T] are highly non-local and not known. Fortu- 
nately, we are still able to classify some phases and phase 
transitions from Eqnj3] without knowing these relations. 
If < 0/ >= for every I = 1,2, 3, the system is in the SF 
state. If < 0; >t^ for at least one I, the system is in 
the insulating state. 

The Harper's equation in the unit cell chosen in the 
Fig. 4 ir*2 was derived in Eqn. A2 hA The authors in 
Refill chose a different unit cell in their Fig.l which is 
the same as Fig.l. In order to be consistent with the for- 
malism developed in=i, we can derive the Harper's equa- 
tion at the filling factor f — 1/q corresponding to the 
unit cell in Fig. lb: 

c a m -!(k x ,k y ) + e tk Hl + e l ^ +2 ^)c a m (k x ,k y ) 

= e(k x , k y )c m (k x , k y ), 

c b m+1 (k x ,k y ) + e-^(l + e-^+ 2 ^)c b m (k x ,k y ) 

= t(k x ,k y )c c ] n (k Xl ky) (4) 

where m — 0, 1, • • • , q — 1 and a, b are two sublattices of 
the honeycomb lattice, —ft/q < k x < n/q, —tt < k y < tt 
are inside the reduced Brillouin zone. The q minima 
of the 2q bands e(k x ,k y ) was found to at (k x ,k y ) = 

(0.27T/0. 
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In the following, we focus on q = 3 case. One ing eigenvectors of the 6x6 matrix at the 3 minima 
only need to find the eigenvalues and the correspond- (k x ,k y ) — (0,27r//3) at / = 0,1,2: 



Ar = 







2pi2-7rZ/3 

l 
o 



i2irl/3 



(1 



i2ir/3\ 



J2ttI/3 



(1 



2e 



-i2irl/3 


1 







-i2nl/3 




1 

(1- 







i2-k/3\ 



(5) 



In fact, we only need the eigenvector c^(Z = 0) = 
c m> c m(' = 0) = c b n at I = 0. The lowest eigen- 
value is e = —(1 + 2cos2-7r/9)£, the corresponding 
eigenvector is (c" n ,cf n ) = [(2cos47r/9 + 2cos27r/9 + 
1), 1/2 - iV3/2,2cos27r/9, (2cos4tt/9 + 2cos2tt/9 + 
l),2cos27r/9, 1/2 - iy/3/2]. Note that c b = c%,c\ = 
c§ , c| — cf . The eigen- vectors at I = 1, 2 can be achieved 
by the magnetic translation^ 2 ! along ci\: c^(Z) = 
c a m Lu- ml ,c b m {l) = c b m u}- ml uj l where w = e l2 * f . The 



eigenfunctions at the three minima (k x ,k y ) = (0, 2nl/3) 
are ^f(x) = E^oC^e^™ 1 ^,^) = 
Em~io c b m (l)e l27Tf ( mai+la2 *> where x = a 1 a 1 + a 2 a2 belongs 
to the sublattice a and x+5, 5 = l/3dl+2/3d2 belongs to 
the sublattice b. One can write the two component vortex 
field at (k x ,k v ) = (0,2^/3) as ^ t (x) = (i/;?(x),ipf(x*)) . 
Then one can write the total vortex field as the expansion 
= Ef=o namely: 



J 



2 

^ a (x) = J2 c„e l2 ™" ai/3 [£o + 6«""e l2 ™ 2/3 + 6w m e" l2 ™ 2/3 ] 

2 

^(f) = ^c^e i27rmai/3 [Co + 6w _m+1 e i27ra2/3 + 6w ro_1 e- i27ra2/3 ] (6) 

m— 



for the vortex fields at sublattice a and sublattice b re- 
spectively. 

In the EqnJHl the £oi£ij£2 are determined by the 
4>q, <j)\, (f>2 in the permutative representation in Eqn.2. 
Plugging in the mean field solutions in Sec. II-B, one can 
determine the vortex fields in the whole dual honeycomb 
lattice. Therefore, one can extract the corresponding va- 
lence bond order in the direct triangular lattice. 

The physical quantities can only be the gauge invariant 
quantities: the densities at different sites in sublattices a 
and b: 

\M*)\ 2 , \Mx)\ 2 (7) 

and the bond quantities between sublattice a and sublat- 
tice b : 

4{x)e a ^ a ^ a {x) = K-iI (8) 

where x belongs to the same unit cell shown in Fig. lb. All 
the other bonds having no phase factor from the gauge 
field. The real part K gives the kinetic energy between 
the two sites. The imaginary part / gives the current 



between the two sites. Note that the "-" sign in Eqn(S] 
is important which make the definition of the current to 
be consistent with that of gauge factors in the Fig. lb. 
X P = J2 P I around any plaquette p ( a hexagon ) gives 
the chirality around the hexagon. The chirality Xp leads 
to a boson density at the center of the hexagon p. 



A. Ising limit v > 0. 

If v > 0, the system is in the Ising limit, only one of 
the 3 vortex fields condense. For example, substituting 
(j>0 = = cj>2 = into Eqns l2l6l7l8l one can evalu- 

ate the vortex densities, kinetic energy and the current 
in the dual honeycomb lattice. For simplicity, we only 
show the currents in the Fig. 2 which are obviously con- 
served. By counting the segments of currents along the 
bonds surrounding the X, Y, Z lattice points, paying spe- 
cial attentions to their counter-clockwise or clock wise 
directions, we can calculate the densities at these points 
n x — 1/3 + 6a;/ > 1/3, n y — n z — 1/3 — 3x1 < 1/3 where 
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(a) (b) 



FIG. 1: (a) Bosons at filling factor / are hopping on a trian- 
gular lattice ( solid line ) which has three sublattices X, Y, Z. 
Its dual lattice is a honeycomb lattice ( dashed line ) which 
has two sublattices a and b. (b) The bond phase factors in 
a dual honeycomb lattice at / = 1/3. The direction of the 
gauge field is important. At f — 2/3, the bond phase factors 
are just complex conjugate of those at / = 1/3. 




FIG. 2: The X-CDW state in the Ising limit v > at / = 1/3. 
At / = 2/3, one can just reverse the current flow and per- 
form a particle and hole transformation in / = 1/3, namely, 
exchange particles ( red dots ) and vacancies ( black empty 
circles ). 



the I is the current flowing around the X lattice site in 
Fig. 2. The x ~ n x — n y > can be thought as a CDW 
order parameter and can be tuned by the distance away 
from the SF to the CDW transition in the Fig. 3a. When 
one tunes the density n y — n z — which stand for vacan- 
cies, then n x — 1 which stands for one boson. This cor- 
responds to the classical limit t = in EqnJT] With the 
quantum fluctuations t > 0, then n x < l,n y = n z > 0. 

Sown in Fig. 2 is a CDW state corresponds to bosons 
occupying one of the three sublattices X, Y, Z in Fig. la, 
so it is 3 fold degenerate. Eqn|3]is an expansion around 
the uniform saddle point < V x A >= / = 1/3 which 
holds in the SF and the VBS ( to be discussed in section 
II-B ) . In the CDW state, the fact that there are currents 
flowing on the dual honeycome lattice indicates that the 
average densities are re-distributed. Because of sharp 
change of the saddle point from EqnJ3] in the SF side 
to Eqn|H]in the CDW side, we expect that the CDW to 
superfluid transition at exactly 1/3 filling in Fig. 3a is first 
order. 

Well inside the CDW phase, there is a large CDW gap 
&-CDW- When studying the physics at slightly away from 
/ = 1/3 filling along the dashed line in Fig. 3a, one has 
to choose a different saddle point where < V x A x >= 
1 — 2a — n x , for the sublattice X and < V x A yz >= a = 
n y = n z for the two sublattices Y and Z should be used. 
The a — > a c < 1/3 limit corresponds to approaching the 



CDW to the SF transition in Fig.3a from the CDW side, 
while the a — > + limit corresponds to the classical limit 
t/V\ — > 0. It is easy to see that there is only one vortex 
minimum (j) yz in such a staggered dual magnetic field 
with a < a c < 1/3, so the effective action well inside the 
CDW state is: 

Ccdw = |(^-^)0 yz | 2 +r|0 yz | 2 + W |^ z | 4 + --- 

+ ife^Mf-W/V) 2 (9) 

where the vortices in the phase winding of <p yz should 
be interpreted as the the boson number. Note that the 
gauge field A x is always massive with a large CDW gap 
AcdWi so was already integrated out in Eqn|9l In the 
direct picture discussed in the online supporting material 
ir>22, the gauge field A v * stands for the uniform density 
fluctuation near q = 0, while the A* stands for the stag- 
gered density fluctuation near Qn — (27r/3, 0) which has 
a large CDW gap in Fig.22. 

Eqn[9] has the structure identical to the conventional 
q = 1 component Ginzburg-Landau model for a " super- 
conductor " in a "magnetic" field. It is easy to see that 
if < cj) yz >y^ 0, the system is in the CDW state where 
both A yz and A x are massive ( Fig.22a). If < 4> yz >= 0, 
the system is in the CDW supersolid ( CDW-SS ) state 
where there is a gapless superfluid mode represented by 
the dual gauge field A v * in Fig. 22b, it still has the gapped 

mode near q — Qn, it also has the same lattice symme- 
try breaking patterns as the CDW state at 1/3 filling. 
Therefore we show that the CDW solid to the CDW- 
SS transition driven by the chemical potential is in the 
same universality class as that from a Mott insulator to 
a superfluid transition, therefore have exact exponents 
z = 1,v = 1/2, rj = with logarithmic corrections ( 
Fig. 3a). It is known the SF is stable against changing 
the chemical potential ( or adding bosons ) in Fig. 3a. 
There must be a transition from the CDW-SS to the SF 
inside the window driven by the quantum fluctuation r 
in the Fig. 3a. The universality class of this transition is 
likely to be first order. 

In the direct lattice picture, we can simply renormalize 
away the sublattice X where the bosons sit and focus on 
the effective boson hopping on the two sublattices Y and 
Z which form a honeycomb lattice whose dual lattice is a 
triangular lattice ( Fig. la). Then the dual vortex action 
in the dual triangular lattice is given by Eqn|9l From 
this direct picture, we can see why a CDW SS can be 
realized much more easily in a triangular lattice than in 
a square lattice: in a triangular lattice, the Y and Z sub- 
lattices still form a connected honeycomb lattices where 
bosons can move without going through the sublattice X 
. While in a square lattice, in a checkboard CDW ( see 
Fig. 19a ), bosons can not hop on sublattice A without 
going through sublattice B. So for a hard core boson, 
checkboard SS is unstable against phase separation be- 
cause bosons can not hop on to the sublattice B. But in a 
soft core case, the checkboard SS could be stabilized be- 
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r<0,v>0 
or r<0,v<0, w<0 

(a) 



CDW-VB-SS < 



r<0, v<0, w>0 



(b) 




FIG. 4: The period 3 stripe phase in the easy plane limit 
v<0,w<0atf — 1/3. At / = 2/3, one can just reverse the 
current flow and perform a particle and hole transformation 
in / = 1/3. 



FIG. 3: Phase diagram of the EBHM in triangular lattice from 
the DVM in this paper, (a) The Ising limit v > 0. There is a 
CDW supersolid slightly above the 1/3 filling which has the 
same lattice symmetry breaking as the C-CDW in the Fig.2. 
The Easy-plane limit v < 0, w < 0. There is a stripe super- 
solid slightly above the 1/3 filling which has the same lattice 
symmetry breaking as the stripe in the Fig. 4. (b) The Easy- 
plane limit v < 0, w > 0. There is a CDW-VB-SS slightly 
above the 1/3 filling which has the same lattice symmetry 
breaking as the CDW-VB in Fig. 7. When the filling factor 
is increased further from the 1/3 filling, there is a first order 
transition from the CDW-SS (VB-SS) to In-commensurate 
CDW ( or In-commensurate VBS ). The IC-CDW or IC-VBS 
are stable only when there are sufficiently long range interac- 
tions. The thin ( thick ) line is the 1st ( 2nd ) order transition. 
As shown in the text, the 1st order transition is a strong 1st 
order one in (a) and a weakly 1st order one in (b). 

cause bosons can hop on the sublattice B too. A striped 
SS in a square lattice could be stabilized even in a hard 
core case, because boson can hop easily along one direc- 
tion, so it is must easier for the bosons to overcome the 
barriers along the other direction to achieve an effective 
hopping. 

B. Easy-plane limit v < 0. 

If v < 0, the system is in the easy-plane limit, all the 
three vortex fields have equal magnitude cf>i — \<j>\e 1 and 
condense. Well inside the VBS side, the mean field solu- 
tion holds. The relative phases can be determined by the 
sign of w. In the easy-plane limit, Eqn|3]is similar to the 
action in Tri-layer quantum Hall systems^,. Introducing 
the center of mass, left and right moving modes as done 
inRefi 2 ?: 9 C = 9 + 9i+6 2> 9i = 6 - 9 2 ,9 r = 9 - 29 x +9 2 , 
EqnJ3] becomes: 

Cvbs = {\d^ c - Atf + \{e^xd v A x -2ir&f5^) 2 + ••• 

30/ . 30; 30 r . , irA 

+ 2w cos — (cos — + cos ) (10) 

In the bipartite lattices in a square and honeycomb 



lattices^, in the easy-plane limit, the system is always 
in a VBS state. As shown in the following, the easy-plane 
limit in a triangular lattice is much involved, in fact, there 
is always a CDW components, so it never corresponds to 
a pure VBS state. 

(1) If w < 0, the mean field solution is 0i — 9q = 
2irm/3,9 2 — 8q = 2im/3 where m,n = 0,1,2. Substi- 
tuting O = 0i = #2 = into Eqns l2l6l7l51 one can 
evaluate the vortex densities, kinetic energy and the 
current in the dual honeycomb lattice. For simplicity, 
we only show the currents in the Fig. 4 which are con- 
served. By counting the segments of currents along the 
bonds surrounding the three rows label, paying special 
attentions to their counter-clockwise or clock wise di- 
rections, we can calculate the densities at these points 
m = 1/3 + Axl > 1/3, n 2 = n 3 = 1/3 - 2x1 < 1/3 
where the / is the current flowing across the lattice. The 
x ~ ni — n 2 > can be thought as a stripe order param- 
eter and can be tuned by the distance away from the SF 
to the stripe transition in the Fig. 3a. When one tunes 
the density n 2 = rc.3 = which stand for vacancies, then 
ni = 1 which stands for one boson. This corresponds 
to the classical limit t = in Eqn[TJ With the quantum 
fluctuations t > 0, then n\ < 1, n 2 — n?, > 0. 

Well inside the stripe phase, it is legitimate to set 
C = 30 o + ^r(m + n)A = -^n,0 r = ^(n - 2m) 
in Eqn llOl the last term reaches its minimum Aw. It 
is this w term which drives the formation of the stripe 
CDW in the Fig. 4. There is a large stripe CDW gap 
^■sCDW ~ w. When studying the physics at slightly away 
from / = 1/3 filling along the dashed line in Fig. 3a, then 
Eqn. [TU] reduces to Eqn|H] Then the discussions follow- 
ing the Eqn|9]also follow with the underlying CDW state 
as the stripe CDW state shown in Fig. 4. The supersolid 
state is a stripe supersolid. The excitation spectra across 
the stripe to the stripe supersolid transition is also given 
by the Fig. 22 with A s cdw ~ w. Alternatively, starting 
from the new saddle point for the dual gauge gauge fields 
corresponding to Fig. 4: < V x A 1 >= 1 — 2a = nx, for 
the row 1 and < V x A 23 >= a = n 2 = 723 for the row 
2 and 3, one can construct the same effective action as 
Eqn|9] after making the replacement x —> l,yz —> 23. 

It is interesting to see that one is not able to get the 
stripe solid state from the DVM in bipartite lattices in 5 ( 
see also appendix B and C ), here one can get the stripe 



FIG. 5: The vortex field of the CDW-VB phase in a tr: 
angular lattice in the easy plane limit v < 0, w > a 
/ = 1/3. Shown is 8 = 0, 0i = 2tt/9, 6» 2 = -2tt/9 case. 
The (0, 0) sets the origin. There are 5 different non-zero mag- 
nitudes shown with 5 different lengths: pi = 3 sin(47r/9), p2 = 
4 sin(7r/3) cos(27r/9) cos(7r/9), p3 = 2sin(7r/3) cos(7r/18), p4 = 
2sin(7r/3) cos(27r/9), ps — sin(7r/3). Note particularly the lat- 
tice points where the vortex fields vanish ! The values of the 
vortex field are listed in the appendix A. 

solid state from the DVM in a triangular lattice from a 
easy plane limit. 

(2) If w > 0, one of the 18 equivalent solutions is 
0- L -6 o = e o -02 = 2tt/9, then 6>i - 2 = 4tt/9. Substitut- 
ing this solution into Eqns l2l6l we find the explicit values 
of the vortex fields at the two sublattices a and b of the 
dual honeycomb lattice graphically shown in Fig. 5 and 
listed in the appendix D. Substituting the vortex field val- 
ues in Eqn IDll into Eqns. !7l8l one can evaluate the vortex 
densities in FigjSl kinetic energy in Fig|B]and the current 
in Fig|7] By counting the number of currents along the 
bonds surrounding the red, green and black lattice points, 
paying special attentions to their counter-clockwise or 
clockwise directions, we can calculate the densities at 
these points n r = 1/3 + 2x[Ii + {I\ — I2)] > 1/3, n g — 
l/3 + 2x[I 2 + (h-h)} > 1/3, rib = l/3-2x[h+h] < 1/3 
where h = sin ^ + sin ?f > 7 2 = sin ~ + sin ?f . The 
x can be tuned by the distance away from the SF to 
the CDW+VB transition in Fig. 3b. When one tunes the 
density n b — 0, then n r — y^pj- > n g — j +j with 
n r /n g = h/h, n r + n g = 1. 

At / = 2/3 case, one only need to change n r —> h r = 
1 - n r = 2/3 - 2x[h + (h - I 2 )] < 2/3, n g -> n g = 
1 - n g = 2/3 - 2x[I 2 + (h - h)} < 2/3, n b -> n b = 
1 - n b = 2/3 + 2x[h + I 2 ] > 2/3. When one tunes 
the density nj, = 1, the n r = j^j^ < n g = 77+77 with 
n r + ri g = 1. When one tunes the density h r — 0, then 

h 9 = 21I-1I >nb = 2w; with n g +n b = 2. 

It is important to observe that as shown in the vortex 
field FigJSJ the vortex fields vanish both at the centers 
of the red loop and the green loop, therefore, both the 
kinetic energy and the current emanating from the cen- 
ters vanish. This fact indicates that there are local SF ( 
or VBS order ) around the these two dual lattice points 
as shown by the red and green triangles in FigJ7] It is 
interesting to compare this with the known fact that in- 
teracting bosons in a kagome lattice at / = l/2^2£ are 
always in a SF state due to the localization of the vor- 



FIG. 6: The Kinetic energy of the CDW-VB phase in a 
triangular lattice in the easy plane limit v < 0, w > 
at / = 1/3. Shown is 8 = 0, 6»i = 2tt/9, 9 2 = -2tt/9 
case. The (0, 0) sets the origin. There are 4 different 
non-zero kinetic energies shown with 4 different colors: 
Kx = 12sin(47r/9)sin(7r/3)cos(27r/9)cos(7r/9), K 2 = 
8sin 2 (7r/3)cos(27r/9)cos(7r/9)cos 2 (7r/18), K 3 = 
4sin 2 (7r/3) cos(2tt/9) cos 2 (tt/18), K 4 = 
2 sin 2 (7r/3) cos(57r/18). Note particularly the bonds where 
the kinetic energies vanish ! These bonds emanate from the 
lattice points where the vortex fields vanish, see the FigfS] 

tices ( a flat vortex band ). The difference is that here 
there is only a local SF around the the centers of the red 
loop and the green loop, while there is a global SF across 
the whole lattice in the latter case. However, the vortex 
field is non- vanishing at the black point, so only the cur- 
rent vanishes, but the kinetic energy as shown in Fig (5] 
does not vanish, it shows there are local CDW around the 
dual lattice point as shown by black dots in the FigtZJ So 
this state has both VBS and CDW which can be called 
CDW-VB state, it is a mixed state unique to a frustrated 
lattice. 

Well inside the CDW-VB phase, it is legitimate to set 
9 C = 30 o ,fy = -^,6 r = -If in EqndUl the last term 
reaches its minimum — |r. It is this w term which drives 
the formation of the CDW-VB in the Fig. 7. There is a 
large CDW-VB gap Aqdw-vb ~ \w\. When study- 
ing the physics at slightly away from / = 1/3 filling 
along the dashed line in Fig. 3b, then Eqn. [TU] reduces 
to Eqn[9] Then the discussions following the Eqn|9] also 
follow with the underlying state as the CDW-SS state 
shown in FigJTJ The supersolid state is a CDW-VB su- 
persolid. The excitation spectra across the CDW-VB to 
the CDW-CB SS transition is also given by the Fig.22 
with Acdw-vb ~ \w\- Alternatively, starting from the 
new saddle point for the dual gauge gauge fields corre- 
sponding to Fig. 7: < V x A r >= n r for the red site, 
< V x A 9 >= n g for the green site and < V x A b >= n b 
for the black site, one can construct a similar effective 
action as EqnO 

In Refi^i, using the density operator formalism on a 
triangular lattice ( See appendix D ), the authors iden- 
tified a bubble phase in Fig. 25 which is a CDW state. 
This bubble CDW phase is completely different from the 
CDW-VB phase in Fig. 7. It is important to identify the 
CDW-VB phase from the the density operator formal- 
ism. Some general problems associated with the density 
operator formalism are examined in appendix D. 



FIG. 7: The CDW-VB phase in a triangular lattice in the 
easy plane limit v < 0, it) > at / = 1/3. Shown is 6a = 
O,0i = 27r/9,6» 2 = -2tt/9 case. The (0,0) sets the origin. 
The red current is I\ — sin ^ + sin , the blue current I2 = 
sin ^ + sin 2^ . The currents are conserved at all the lattice 
points. At / = 2/3, one can just reverse the current flow and 
perform a particle and hole transformation in / = 1/3. For 
6*o = O,0i = — 27r/9,02 = 27r/9, one only need to rotate the 
figure by 180°. 



We conclude that the transitions from the X-CDW, 
stripe-CDW, CDW-VB solid to the corresponding super- 
solids driven by the chemical potential in the Fig. 3 is also 
in the same universality class as that from a Mott insu- 
lator to a superfmid. 



C. Valence bond state and absence of Valence 
bond supersolid in a triangular lattice 

It is easy to draw the VBS state in a triangular lat- 
tice in the Fig. 8 This is a plaquette state similar to the 
plaquette state in a square lattice shown in Fig. 20b In 
contrast to the dimer VBS in a square lattice, there is no 
dimer VBS at a triangular lattice. So the uniform saddle 
point < V x A >= f = 1/3 holds in both the SF and 
the VBS. As shown in the previous two sections, from 
the dual vortex theory, one is not able to get a pure VBS 
state at / = 1/3 shown in Fig. 8 in either the Ising limit 
or easy-plane limit. Fig. 7 contains the VB component, 
but also a CDW component. Because the excitation of 
this state necessarily involves a density excitation shown 
in Fig. lib, drawing the insights gained from the direct 
first order transition from a pure VB state to a SF in 
a Kagome lattice Fig. 12 to be discussed in Sec.III-2, we 
conclude that there is no VB-SS in a triangular lattice, 
in sharp contrast to bipartite lattices discussed in* and 
in the appendix B and C, there can only be a direct first 
-order transition from the VB to the SF. 



D. Implications on QMC simulations in triangular 
lattice 

(1) Nearest Neighbor (NN) EBHM in a triangular lat- 
tice 

The simplest EBHM in EqnfT] is the hard core bosons 
hopping in a triangular lattice with only nearest neighbor 



(a) (b) 

FIG. 8: Two kinds of valence bond solids (VBS) in a triangle 
lattice. ( a ) Uniform VBS state, (b) Staggered VBS state. 

(nn) interaction U = 00, Vi > 0: 

H nn = -t ^ (b}bj+h.c.)+Vi ^ riiUj - fi^^m (11) 

<ij> <ij> i 

This model was studied expensively by QMC ii>L2r— 
where the hard core constraint rii = b\bi < 1 consider- 
ably reduces the size of Hilbert space in QMC simula- 
tions. The QMC results ir*i2i2& leads to Fig. 9a. When 
comparing the Fig. 3a In the Ising limit r < 0, v > 
achieved in this paper with the Fig. 9a by QMC, one can 
see that the solid state at 1 /3 is just the X- solid state in 
Fig.2. the SS-i is just the CDW-SS. The IC-CDW state in 
Fig. 3a can not be realized in Fig. 9a, because it can only 
be stabilized by much longer-range interactions than V\ . 
The dual vortex effective action to describe the transi- 
tion from the SF to the 6-fold solid is given by EqnJ3] 
in the Ising limit r < 0, v > 0, it was found to be very 
weak first order by QMC i n 19 i 20 . The effective action to 
describe the transition from the 6-fold solid to the stripe 
supersolid slightly away from 1 /3 is given by Eqnj9l How- 
ever, the universality class of the 2nd phase transitions 
in Fig. 2a was not studied i n 19 i 20 . Comparing the QMC 
results with Fig. 3a from the DVM, we can see the solid 
at 1/3 filling to the interstitial induced supersolid ( SS-i 
) slightly above 1/3 filling is indeed second order transi- 
tion in the same universality class as that from a Mott 
insulator to a superfluid. However, the SF to the SS-i 
is likely to be weakly first order. Note that sometimes 
it is difficult to distinguish a 2nd order transition from 
a weakly first order one from QMC. More refined QMC 
simulations are needed to address the nature of the SF 
to SS-i transition. 

(2) Next Nearest Neighbor (NNN) EBHM in a trian- 
gular lattice: Stripe phase, Stripe-SS phase 

The next simplest EBHM in EqnfT] is the hard core 
bosons hopping in a triangular lattice with NN and next 
nearest neighbor (nnn) interactions U = 00, Vi > 0, V% > 
0: 

H nnn = -t 2J {b\bj + h.c.) + Vi 2J n % n i 
<ij> <ij> 

E n » ( 12 ) 

«ik» i 
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It was studied by QMC in22. A period-3 striped soild 
state is found at / = 1/3 and a stripe-supersolid was 
found at slightly large than 1/3. The QMC results irM 
leads to Fig. 9b. When comparing the Fig. 3a in the easy 
plane limit r < 0,v < 0,w < achieved in this paper 
with the Fig. 9b by QMC, one can see that the stripe 
state at 1/3 is just the stripe state in Fig. 4. The dual 
vortex effective action to describe the transition from the 
SF to the stripe solid is given by Eqnj3] in the easy plane 
limit r < 0, v < 0, w < 0, it was found to be strongly first 
order by QMC in 23 . The effective action to describe the 
transition from the stripe solid to the stripe supersolid 
slightly away from 1/3 is given by Eqn[9] However, the 
universality class of the 2nd phase transitions in Fig. 9b 
was not studied either ir* 2 ^. 

(3) Searching for the CDW-VB phase in the NNN 
EBHM model in a triangular lattice 

A meta-stable bubble solid phase in Fig. 25 was also 
found in the QMC ir* 2 ^. This bubble solid phase has 
higher energy than the stripe solid phase. The very in- 
teresting CDW-SS phase in Fig. 7 was not searched in the 
QMC in^ in any parameter regimes. In order to identify 
this phase, in addition to the QMC calculations of the 
superfluid density and the density structure factor in 2 ^, 
a bond structure factor— need also be studied to search 
for the VB ordering in the Fig. 7. 



III. SOLIDS AND SUPERSOLIDS IN KAGOME 
LATTICE NEAR 1/3(2/3) 

The Kagome lattice ( Fig. 10a ) has 3 sublattices a, b 
and c, its dual lattice is the dice lattice ( Fig. 10a). Due 
to too strong quantum fluctuations, as shown in 2 - the 
spin wave expansion does not work anymore in Kagome 
lattice. Here I am trying to investigate the phases and 
quantum phase transitions when slightly away from the 
1/3 filling from the DVM. Because of the P-H symmetry 
at U — oo, the results are equally valid near 2/3 filling. 
As shown in^ at q = 2, because the lowest dual vortex 
band in the dice lattice is completely flat, so the dual vor- 
tices are completely localized, so it has to be a superfluid, 
in sharp contrast to q = 2 on triangular lattice Fig. 9. For 
q = 3, from Table 1 of^, we can see there are two minima 
at (0,0) and (— 27r/3, 27r/3). Let's label the two eigen- 
modes at the two minima as — 0,1. The effective 
action invariant under all the MSG transformations upto 
sixth order terms was written down in 2 ^. Again, inside 
the superfluid phase, moving slightly away from the 1/3 
filling / = 1/3 corresponds to adding a small mean dual 
magnetic field Sf = f — 1/3 in the action derived in2&. 
Upto sixth order, the action is Lsf = Cq + C\ + C 2 '- 

£0 - ^|(a M -a M )0 ; | 2 + ^l^| 2 + ^(e^A^A A -27r^/,5 Al 
I 

C 2 = wUl^f + h.c] ( 




(a) (b) 

FIG. 9: (a) QMC results ini&2£ on triangular lattice with 
U = 00, Vi > 0. Slightly above ( below ) the 1/3(2/3) filling, 
there is a SS-i ( SS-v) which has the same lattice symmetry 
breaking patterns as the solid at 1/3 ( 2/3 ) where the bosons 
occupying ( un-occupying ) one of the 3 sublattices X, Y, Z in 
Fig. la. SS-i and SS-v are related by the P-H transformation. 
Exactly at half filling where /i/Vi = 3, SS-i and SS-v coexist 
at any possible ratio, there is a first order transition from 
SS-i to SS-v driven by the chemical potential [i across the 
thick dashed line, (b) Suggested phase diagram on triangular 
lattice with U = 00, Vj = 0, V 2 > by the DVM in this 
paper. The Solid at 1/3 is the period-3 stripe phase studied 
numerically first The SS-i slightly above 1/3 filling is 

the corresponding period-3 stripe supersolid phase studied by 
recent QMC in^. If the solid at 1/3 is the CDW-VB solid 
phase in Fig. 7, then the corresponding possible CDW-VB SS 
phase maybe unstable against phase separation. The thin ( 
thick ) line is the 1st ( 2nd ) order transition. As explained 
in the text, the first order transition is a strongly first order 
one. The universality class of the 2nd order phase transitions 
was not studied 



where is a non-compact U(l) gauge field. 

In fact, the form of this action is exactly the same as 
that in honeycomb lattice at q — 2 first derived ir£. 
However, the physical meanings of the two eigenmodes 
4>i,l = 0,1 are completely different in the two lattices, 
this fact leads to completely different physics in the two 
lattices. 

Ref.— studied the vortex environments in the dual dice 
lattices, but did not evaluate any gauge invariant physical 
quantities. The vortex enurements are gauge dependent, 
so may not be used to characterize the symmetry break- 
ing patterns. Ref.— studied the bosons hopping on a Dice 
lattice at an integer filling. These bosons have only on- 
site interaction, but are subject to some artificial gauge 
fields which can be generated in cold atom experiments. 
They identified the vortex fields and the gauge invari- 
mit currents at / = 1/3. These artificial gauge fields are 
static. While the gauge field in all the dual vortex ac- 
tions Eqn ll3ll6ll4l on the dual Dice lattice are quantum 
fluctuating which are responsible for the very existences 
3>T the CDW-VB-SS. 
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(a) 0>) 



FIG. 10: ( The VBS around the up triangle ( triangular va- 
lence bond (TVB) ) in a Kagome lattice in the Ising limit 
v > at / = 1/3. The thick triangle means one at / = 1/3 
( or two at / = 2/3 ) boson(s) is (are) moving around the 
up triangle, (b) The saddle point of the dual gauge field in a 
magnetic field unit cell in the dice lattice. The saddle point 
stays the same as that in the SF phase. 



A. Ising limit, v > 0. 

If v > 0, the system is in the Ising limit, only one of the 
2 vortex fields condense. Let us take (0o) 0> (<f>i) = 0- 
Using Eqns l7l8l one can evaluate both the vortex density, 
the kinetic energy and current in the dual dice lattice. 
This has been done in2£ in a very different context. There 
is no current flowing in the dual lattice which shows there 
is no CDW order. All the densities at the direct lattice 
are fixed at 1/3. It is important to point out that the 
vortex fields at the dual center of the red triangles in 
the Fig. 10a vanish, so both currents and the kinetic en- 
ergies emanating from these dual centers are vanishing, 
this fact indicates that there is a local VBS order around 
red triangles as shown in the Fig. 10. This is just oppo- 
site to bipartite and triangular lattices discussed in£ and 
the last section where the Ising limit corresponds to the 
CDW cases. This crucial difference is responsible for the 
absence of supersolid in the Ising limit as shown in the 
following. Because the saddle point structure stays the 
same across the SF to the triangle- VBS ( TVB ) transi- 
tion at / = 1/3 as shown i the Figl2a, Eqn|T3] still holds 
in the TVB side, so the transition is a weak first order 
one which breaks both the U(l) and the Z-i exchange 
symmetry between cj>g and tf>\. In contrast, in the Ising 
limit of both the bipartite and triangular lattice, different 
saddle points for the dual gauge fields need to be chosen. 

Slightly away from 1/3 filling, there must be a direct 
first order transition from the TVB to the superfiuid as 
shown in Fig. 12, there is no TVB supersolid intervening 
between the RVB and the SF, in sharp contrast to the 
cases in bipartite and triangular lattices. 

Note that due to the lack of dimer VBS on the Kagome 
lattice at / = 1/3, any excitation above the TVB state 
in the Fig. 10 will necessarily involve the breaking the the 
TVB into a density excitation shown in Fig.ll. This gap 
is of density origin represented by the gauge field fluctu- 
ation in Eqn ll3l This is in sharp contrast to the dimer 
VB or plaquette VB state in bipartite lattices such as 
square and honey come lattice. In these bipartite lat- 
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FIG. 11: A local density excitation above the TVB state 
in Fig. 10. Its translational moving through the whole lattice 
leads to the excitation spectrum in the upper right Fig. 12. 
Compare with the local excitation in a VBS in Fig. 19. 
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Ising limit 

FIG. 12: The phase diagram driven by the chemical potential 
fi. There is a direct first order transition from the TVB to 
the superfiuid. There is no supersolid intervening between 
the TVB and the SF in any case. Upper right is the gapless 
excitation due to the gauge field A M in the SF. Lower right 
is the gapped excitation due to the Higgs mechanism of the 
gauge field A M inside the TVB. 

tices, there are VBS excitation gaps due to the operator 
A cos which only involves dimer flips instead of any 
density fluctuation shown in Fig. 21. Here, such an oper- 
ator is lacking in the Ising limit, so there is no such VBS 
gap. 

B. Easy-plane limit v < 0. 

If v < 0, the system is in the easy-plane limit, Eqn|3]is 
similar to the action in Bi-layer quantum Hall systems^, 
the two vortex fields have equal magnitude <j>\ = \4>\e 101 
and condense. The relative phases can be determined by 
the sign of vA. When moving from 1/3 to 1/2 where the 
bosons has to be in a SF state, Eqn[l6]can be rewritten 
as: 

Cvbs = (~d»9 + - \) 2 + \(^xd,A x -2irSf5^) 2 + --- 
+ \i d ^-f + 2wcos36C (14) 
where 8± = 0q ± 0\ . 
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FIG. 13: The charge density wave + VBS state in a Kagome 
lattice in the easy-plane limit v < 0, w > at (a) / = 1/3 
(b) / = 2/3. The red solid circle is a boson, the empty red 
circle is an vacancy, (a) and (b) are related by time reversal 
symmetry and particle-hole symmetry. 



In the following, we discuss both cases separately. 
(l)w>0. CDW+VBS order: 

When w > 0, 6L = 9 - B x = (2n + l)vr/3, using 
Eqns l7|8l one can evaluate both the vortex density, the 
kinetic energy and current in the dual dice lattice. This 
has been done ir^ in a different context. The current was 
shown in the Fig ll3l It is important to point out that 
the vortex at the dual center of the red hexagon vanishes 
which indicates the there is a local VB order around the 
dual center as shown by the red hexagon in the FigfT3l 
Again, it is interesting to compare this with the known 
fact that interacting bosons in a kagome lattice at / = 
1/2^2^ are always in a SF state due to the localization of 
the vortices ( a flat vortex band ). The difference is that 
here there is only a local SF around the the center of the 
red hexagon, while there is a global SF across the whole 
lattice in the latter case. By counting the segments of 
currents along the bonds surrounding the red and black 
lattice points, paying special attentions to their counter- 
clockwise or clock wise directions, we can calculate the 
densities at these points n& = 1/3 — 4x1 < 1/3, n r = 
1/3 + 2x1 > 1/3 at the filling factor / = 1/3 where 
the I is the current flowing around the black lattice site 
in FigUS] The x ~ n r — n& > can be thought as 
a CDW-VB order parameter and can be tuned by the 
distance away from the SF to the CDW-VB transition in 
the Fig. 3b. When one tunes the density tt.^ = which 
stand for vacancies, then n r = 1/2 which stands for 3 
bosons are hopping around the the red hexagon to form 
a local VB order. Very similar calculations can be made 
for the filling factor f — 2/3 which is related to / = 1/3 
by a particle-hole transformation shown in Fig. 13b. 113b . 

Well inside the CDW-VB phase, it is legitimate to set 
6_ = 6 - 6% = (2n + 1)tt/3 in EqniTil the last term 
reaches its minimum — 2w. It is this w term which drives 
the formation of the CDW-VB phase in the Fig. 10. There 
is a CDW-VB gap Acdw-vb ~ W. When studying the 
physics at slightly away from / = 1/3 filling along the 
dashed line in Fig. 3a, then Eqn. [13] reduces to EqnfMl 

Alternatively, one can start from the new saddle point 
for the dual gauge gauge fields corresponding to Fig. 14b 
which is < V x A C dw >= n b = 1 - a, < V x A V bs >= 



n r = (1 + a)/2 as shown in Fig, 3b. The lowest energy 
dual vortex band in such a saddle point is found to be: 



E(k u k 2 ) 
A(ki,k 2 ) 



-V64 
cosfci - 
cos(fci 



2A{k!,k 2 ), 
cos(fci — an) 
- k 2 ) — cos(fci 



- 7T < fel, k 2 < 7T 

■ cos k 2 — cos{k 2 + an) 
k 2 - 2an) (15) 



It is easy to see that there is a symmetry E(ki,k 2 ) — 
E(k 2 + n(l + a),ki + n(l-a)). For < a < 0.404, there 
are two minima at k-i — an 12 — arcsin( V. n ^ a 7^ 2 ) ) , k 2 — 
— fci and fci = a7r/2 + arcsin( s ™^ a ^ 2 ) ) + tt, fc? = — fci. 

1 1 i v 2 cos(a7r) ; ' z 1 

The two minima are related by the symmetry. 
When a = 1/3 < 0.404, it reduces to q — 



where there are two minima at (0, 0), (— 



3 case 
which 
> 1/3" 



is the two minima at the superfluid side. The a 
case corresponds to approaching to the transition from 
the CDW+VB to the SF in Fig.3b from the left hand 
side. For a — which corresponds to t/V\ limit in 
Fig. 17a, A(ki,k 2 ) — 2cos(fci — k 2 ), then the minima is 
along the line k\ — k 2 , the band is flat along this direc- 
tion. At q = 2, the band is completely flat as shown in£. 
In fact, when a = 0, all the 6 parallelograms around the 
center in the Fig. 14b have q = 2, the other 6 have q = 1, 
so it is not a surprising that the band is flat along one 
direction. Of course, for a = 1, it reduces to the zero 
magnetic case q = 1 where there is only one minimum at 
(0,0). 

Now for < a < 1/3 , because the two min- 
ima are smoothly connected to those in the SF side, 
let's still label the two minima by <fii,l = 0,1. Then 
the MSG inside the SF state T^,T 2 ,R n / s ,I\ is reduced 
to Ts 1+ s 2 , i?^/3, Iax+a-2 inside the CDW+VBS state in 
Fig. 14. Note that only the translational symmetries 
are reduced, but the rotational symmetry remains the 
same as inside the SF state. The effective action slightly 
away from 1/3 invariant under the MSG transformations 
inside the CDW+VBS state upto sixth order terms is 
^-cdw+vbs = Cq + C\+ C 2 : 



Co = £>P 



VBS\ 



c 2 



"(I0ol 



w 



»o<Pi 



-%A 

f I01l 2 ) 2 - 
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( c a A VBS 



2716/6, 



2\2 



m \ 2 -\m 



where A^ BS is a non-compact U(l) gauge field. 

Note that Eqn[T3] in the SF side and EanJIrjl in the 
CDW+VB side take the same form despite the reduced 
translational symmetry in the CDW+VB state. Note 
that A^ DW is always massive, so was already integrated 
out in Eqn llGI Because the system is already in the 
CDW+VBS, so we expect f < and v < 0. Then we 
can set <fii = |0|e l8i . The relative phases can be deter- 
mined by the sign of w which has the same sign as w > 0. 
Because of the very weak relevance of w term in EqnsfT3l 
and [TBI so we expect the SF to the CDW+VB transition 
driven by r is very weakly first order. This is in sharp 
contrast to the strong 1st order CDW to the SF tran- 
sition in a triangular lattice in Fig. 3a and Fig. 9a. This 
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FIG. 14: ( a ) The CDW+VBS in a Kagome lattice at / = 2/3 
identified in^. The dot means one ( or I— a ) boson occupying 
the site. The O inside the hexagon means 3 ( or 3(1 + a) 
) bosons hopping around the hexagon. It has 3 sublattices 
labeled as a,b,c. ( b ) The saddle point of the dual gauge 
field in a magnetic field unit cell in the dice lattice, a = 
corresponds to (a). The Si and 3,2 are the two basis vectors 
in the dual dice lattice 



is exactly what were observed in the QMC in triangular 
lattic o 19 ' 20 and in kagome lattice^. 

When moving from 1/3 to 1/2 where the bosons have 
to be in a SF state, EanfTfUalso reduces to EqnfMl The 
phase diagram may be similar to Fig. 3b in the triangu- 
lar lattice. The SS slightly away from 1/3 filling has the 
same CDW and VB order as the solid at 1/3, so let's call 
this kind of novel SS as CDW-VB-SS. Then the transi- 
tion from the CDW+VB to the CDW-VB-SS transition 
driven by the chemical potential can be similarly dis- 
cussed as the VBS to the VB-SS in bipartite lattices^, 
so it is is the same universality class as that from the 
Mott insulator to the SF. The CDW-VB-SS to the SF 
transition is 1st order. Although the CDW-VB-SS is not 
stable in the hard core case, it is stable in the soft core 
case. Of course, there is no P-H symmetry anymore in 
the soft core case. This study may inspire more accurate 
QMC to search for this novel CDV- VB-SS state, espe- 
cially in the soft-core case. The excitation spectra across 
the CDW-VB to the CDW-VB supersolid transition is 
also given by the Fig. 22 with Acdw ~nor equivalently 
by Fig. 25 with Avbs ~ w. 

However, Eqn ll6l may break down as a — > which 
corresponds to t/Vi — > limit, because the two minima 
become very shallow. At a — 0, the band is flat along 
k\ = k2, so we expect the regime of the CDW + VB-SS 
shrinks to zero as t/V\ — > as shown in Fig. 17a. 

(2) w < case: 6- fold CDW order: 

When w < 0, 0_ = 6 - B x = 2nn/3, using EanslTISl 
one can evaluate both the vortex density, the kinetic en- 
ergy and current in the dual dice lattice. This has been 
done ir*^ in a very different context. The current was 
shown in the Fig |15l It is important to point out that 
the vortex at the dual center of the black spots is non- 
vanishing, so only the current is zero due to the can- 
celation of the current flowing on the neighboring black 
points, but the kinetic energy is positive. This fact in- 
dicates the there is a local CDW order around the dual 
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FIG. 15: The 6-fold CDW state in a Kagome lattice in the 
easy-plane limit v < 0, w < (a) / = 1/3 (b) / = 2/3. 
The red solid circle is a boson, the empty red circle is an 
vacancy, (a) and (b) are related by time reversal symmetry 
and particle-hole symmetry. 

center as shown in the Fig. 1151 By counting the segments 
of currents along the bonds surrounding the red and black 
lattice points, paying special attentions to their counter- 
clockwise or clock wise directions, we can calculate the 
densities at these points n r = 1/3 + Axl > 1/3, n& = 
1/3 — 2x1 < 1/3 where the I is the current flowing around 
the red lattice site in FigfT51 The x ~ n r — ri\, > can 
be thought as a CDW order parameter and can be tuned 
by the distance away from the SF to the CDW transi- 
tion in the Fig. 3a. When one tunes the density n.& = 
which stand for vacancies, then n r — 1 which stands for 
one boson. This corresponds to the classical limit t = 
in EqnJT] With the quantum fluctuations t > 0, then 
n r < l,n b > 0. In fact, the 6-fold CDW state Fig.15 
has the same symmetry as the CDW-VB state in Fig. 13. 
The only difference is the local CDW order versus a local 
VBS order. 

Well inside the 6-fold CDW phase, it is legitimate to 
set 6*_ = 8q — 8i = 2mr/3 in EqnfT^l the last term reaches 
its minimum 2w. It is this w term which drives the for- 
mation of the 6-fold CDW in the Fig.15. There is a 6-fold 
CDW gap Acdw ~ w. When studying the physics at 
slightly away from f — 1/3 filling along the dashed line 
in Fig. 3a, then Eqn. Q2] with w < reduces to EanlTH 
with w < 0. The discussions following the EqnQj] also 
hold with the underlying CDW state as the 6 fold CDW 
state shown in Fig.15. The supersolid state is a 6 fold 
CDW supersolid. The excitation spectra across the 6- 
fold CDW to the 6-fold CDW supersolid transition is 
also given by the Fig. 22 with Acdw ~ w. Alternatively, 
starting from the new saddle point for the dual gauge 
gauge fields which is identical to that of the CDW-VB 
phase shown in Fig. 14b, one can also derive Eqn ll4l and 
use its following discussions. 



C. Stripe-CDW order: 

As shown in the previous two sections, one is not able 
to get the stripe state at / = 1/3 in Fig. 16 from the 
dual vortex theory in either the Ising limit or easy-plane 
limit. But it is very easy to get this state from the NNN 
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FIG. 16: ( a ) The stripe CDW in a Kagome lattice at / = 1/3. 
The dot means one ( or a > 1/3 ) boson occupying a site in 
lattice a. ( b ) The saddle point of the dual gauge field in 
a magnetic field unit cell in the dice lattice. Compare the 
saddle point with Fig. 14b 

EBHM Eqn|T2]in a Kagome lattice. If we add very strong 
nnn interaction Vi , the bosons will go to the nnnn which 
belong to a again ( Fig. 16), then the bosons will simply 
take one of the 3 sublattices such as a to form a CDW 
alone ( Fig. 16). Slightly away from 1/3, there is a stripe 
SS as shown in Fig. 17b. 

In the direct lattice Fig. 16a, even if taking away the 
sublattice a, the bosons can still move easily along the 
chain bcbc..., so they need only overcome the barrier along 
the chain baba... to achieve an effective hopping along 
baba.... So the superfluid stiffness along the be chain is 
larger than that along the ba chain. In the stripe CDW 
state in Fig. 16a, a different saddle point where < V x 
A a >= l-2a for the sublattice a and < Vx^ >= a for 
the two sublattices b and c should be used. The a — > a c < 
1/3 limit corresponds to approaching the stripe CDW 
to the SF transition in Fig. 17b from the stripe CDW 
side, while the a — > + limit corresponds to t/Vi —> 
limit in Fig. 17b. It is easy to see that there is only one 
vortex minimum fac hi such a staggered dual magnetic 
field with a < a c < 1/3, so the effective action inside 
the stripe CDW state should be the same as Eqn. 
with yz — > be. Because of sharp change of the saddle 
point from EqnJT3]in the SF side to Eqn|4]in the stripe 
CDW side, so the transition from the SF to the stripe 
CDW along the horizontal axis in Fig. 17b is likely to be 
a strong first order. Then the stripe-CDW to the stripe- 
SS transition driven by the chemical potential along the 
vertical axis is s second order transition 



D. Implications on QMC simulations in Kagome 
lattice 

(1) Nearest Neighbor (NN) EBHM in a Kagome lattice 
Eqn[T]in Kagome lattice with U = oo, V\ > 0, namely, 
EqnJTT] in a Kagome lattice, was studied by two recent 
QMC simulations 2 ^ 2 !. It was found that the solid state 
at exactly 1/3 filling has both CDW and VB order. 
The VB order corresponds to boson hopping around a 
hexagon ( shown in Fig. 13b ). This is a unique feature 
of the Kagome lattices. This feature can be intuitively 



understood without any calculation: with very large V\ 
interaction at 1/3 filling, if a boson taking sublattice c, 
then all its 4 nn which belong to sublattices a and b are 
excluded, so bosons will go to the 4 nnn which still be- 
long to a and b instead of c ( Fig. 16a). This is in sharp 
contrast to triangular lattice at 1/3 filling in Fig. la. So 
bosons may move around the hexagon to form some VB 
order in addition to the CDW order. This state is iden- 
tical to the CDW-VB state in Fig. 13a. Although the 
QMC in Ref. 27 found the the SF to the CDW+VB tran- 
sition is a 2nd order transition through a noval decon- 
fined quantum critical point 2 !, the QMC in Ref. 2 '' found 
a weakly 1st order transition from a double-peaked his- 
togram of the boson kinetic energy. The dual vortex 
effective action to describe the transition from the SF to 
the CDW+VB is given by Eqn[T2]in the easy plane limit 
r < 0,v < 0,w > 0. As explained below EqnJTBl it should 
be a weakly 1st order transition from the CDW-VB solid 
to the SF which is consistent with the conclusion in 2 !. 

The effective action to describe the transition from the 
stripe solid to the stripe supersolid slightly away from 1 /3 
is given by Eqn ll6l Drawing the insights gained from 
the square lattice! , we conclude that the CDW-VB- SS is 
unstable against the phase separation in the hard core 
limit studied i n 26 i 27 . But it will be stable in the soft core 
case. 

(2) Next Nearest Neighbor (NN) EBHM in a Kagome 
lattice 

It may be interesting to do QMC with U = oo, V\ > 
0, V<z > to test Fig. 17b. The hamiltonian is the same 
as EqnfT2l but in a Kagome lattice. As elucidated in the 
last paragraph, when both V\ and Vi are large, a CDW 
state where bosons simply take one of the 3 sublattices 
can be stabilized. Because the nn and the nnn neighbors 
are similar as can be seen from Fig. 6a, the vertical axis in 
Fig. 17b need to be replaced by \ij (V1 + V2), the horizontal 
axis stay as t/Vi at fixed large V% ( Fig. 17b ). While 
the conclusion achieved in! that the system has to be a 
superfluid state at 1/2 remains robust. So we conclude 
that when moving from 1/3 to 1/2, the system goes CDW 
to stripe-SS to the SF transition. The first is the second 
order, the second is a first order transition ( Fig. 17b). 
Again, longer-range interactions favor the description in 
terms of the dual vortices and the stability of SS, so the 
Stripe-SS will surely be stable even in the hard core case 
( Fig.l7b ). 

(3) Searching for the TVB state: possible ring ex- 
change interactions 

So far, the TVB state has not been identified in 
any QMC simulations. The ring exchange interaction 
— Kk J2ijki {Kbjbjfii + h.c.) where i, j, k, I label the 4 cor- 
ners of a bow tie in the Kagome lattice is needed maybe 
stabilize this TVB phase. 
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FIG. 17: ( a ) Suggested zero temperature phase diagram in 
Kagome lattice U = oo, Vi > 0. In the soft core case U < oo, 
there maybe a narrow window of CDW-VB-SS sandwiched 
between the CDW-VB solid and the SF. Of course, there is 
no P-H symmetry anymore in the soft core case. In the hard 
core case studied i n 26 ' 27 , the CDW-VB-SS is unstable, there is 
a direct weakly first order transition from the CDW-VB solid 
to the SF. Thin (thick) line is 2nd (first ) order transition. It 
was shown in 3 that at half filling / = 1/2 where n/V\ = 2, 
it can only be a SF state, in sharp contrast to q — 2 on 
triangular lattice Fig. 9. (b) Very similar phase diagram also 
works for U = oo, Vi > 0, Vi > 0. In this case, Vi in (a) need 
to be replaced by Vi + V2, the solid at 1/3(2/3) is simply a 
CDW where the bosons occupying ( un-occupying ) one of the 
3 sublattices a,b,c in Fig. 16a. The SS is simply a stripe SS 
where the superfluid stiffness p s \ 2> psi- The Stripe-SS will 
surely be stable even in the hard core case. As shown in the 
text, the first order transition in (a) is a weakly first order 
one, but strongly first order one in (b) 



IV. CONCLUSIONS 

In summary, we used the DVM developed infill 5 , to 
study quantum phases, especially various kinds of insu- 
lating phases and phase transitions in the two most com- 
mon frustrated lattices such as triangular and kagome 
lattices at and slightly away 1/3 and 2/3 fillings. At the 
commensurate fillings 1/3(2/3), we developed a system- 
atic way to identify the symmetry breaking insulating 
states uniquely and completely by using gauge invariant 
vortex density, kinetic energy and current on the dual 
lattices. We reproduced several known phases found by 
the previous methods, also identified several new phases. 
We contrast different symmetry breaking patterns of the 
insulating phases in both Ising and Easy plane limits in 
the two lattices. By analyzing carefully the saddle point 
structures of the dual gauge fields in the translational 
symmetry breaking side and pushing the effective action 
slightly away from the commensurate fillings, We classi- 
fied all the possible types of supersolids and analyze their 
stability conditions. By using these effective actions, we 
studied different kinds of solids, different kinds of su- 
persolids and the universality class of quantum phase 
transitions from the solids to the supersolids. We also 



compared our results achieved by the DVM with some 
available QMC results and make implications on pos- 
sible future QMC simulations. The results achieved in 
this paper should have direct impacts on ultra-cold atoms 
loaded on optical lattices. Although supersolids on lat- 
tices are different than that in continuous systems such as 
4 He~ or exciton supersolid in electron-hole bilayer sys- 
tem ir>2£, they may share some common properties and 
shed considerable lights on each other. 

The DVM is a magnetic space group ( MSG )i 
symmetry-based approach which, in principle, can be 
used to classify all the possible phases and phase transi- 
tions. But the question if a particular phase will be stable 
or not as a ground state depends on the specific values 
of all the possible parameters in the EBHM in Eqn|TJ 
so it can only be addressed by a microscopic approach 
such as Quantum Monte-Carlo (QMC) simulations. The 
DVM can guide the QMC to search for particular phases 
and phase transitions in a specific model. Finite size 
scalings in QMC in a microscopic model can be used to 
confirm the phases and the universality classes of phase 
transitions discovered by the DVM. The two methods 
are complementary to each other and both are needed to 
completely understand phases and phase transitions in 
EqnH] 

In the square lattices, so far, there are several estab- 
lished examples of nice comparisons between the phe- 
nomenological DVM and the QMC simulations on a mi- 
croscopic model: (1) The transition from the SF to the 
X-CDW ( Fig. 19a ) in the Ising limit v > in a square 
lattice was studied by the DVM- and the QMC on the 
EBHM Eqn[Tl]with NN interaction in a square lattice at 
/ = 1/2 with both hard 3 ^ and soft core^ 5 . bosons. (2) 
The transition from the SF to the X-CDW ( Fig. 19b ) 
in the Ising limit v > in a honeycomb lattice was stud- 
ied by the DVM 5 and the QMC on the EBHM Eqn[TT] 
with NN interaction in a honeycomb lattice at / = 1/2 
with both hard and soft core bosons^. Unfortunately, 
there are still many cases the DVM and QMC can not 
be compared with each other. For example, the DVM is 
not able to get the stripe-solid phase studied extensively 
by QMC 34 on the EBHM of hard core bosons with NNN 
interaction EqnQjJat / = 1/2. Interestingly, the stripe 
solid phase Fig. 4 in a triangular lattice has been achieved 
by the DVM in the easy plane limit v < 0, w < 0. How- 
ever, the stripe solid phase Fig. 16a can not be achieved 
from the the DVM either. It remains an open problem 
how to achieve the strip solid phases in a square and 
kagome lattices by the DVM. On the other hand, so far, 
it has not been possible to study all the VBS phases by 
QMC in some microscopic models in the easy-plane limit 
of the DVM listed in the Fig.20 and Fig.23. The above 
two connections are established only in the Ising limit. 
This maybe due to the fact that a ring exchange inter- 
action may be needed to stabilize all these VBS phases. 
It is still not known what are the specific forms of these 
ring exchange interactions. 

In the frustrated lattices, so far, there are also sev- 
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eral established examples of nice comparisons between 
the phenomenological DVM and the QMC simulations 
on a microscopic model. For example, The DVM 3 -^ pre- 
dicted that the Kagome lattice at 1/2 must be a super- 
fluid due to the complete flat band of the dual vortices. 
No insulating phases can be found at 1/2. This has been 
confirmed by QMC2&. In the following, we list the con- 
nections directly relevant to this paper: (1) The transi- 
tion from the SF to the X-CDW ( Fig. 2 ) in a triangular 
lattice at / — 1/3 was studied by the DVM in^i in the 
Ising limit v > of the effective action Eqn(3]and also by 
the QMGii^ on the EBHM of hard core bosons with NN 
interaction EqnfTTl (2) The transition from the SF to the 
Stripe-CDW ( Fig. 4) in a triangular lattice at / = 1/3 
was studied by the DVM in^i in the easy plane limit 
v < 0, w < of the effective action Eqn|3] and also by 
the QMG^on the EBHM of hard core bosons with NNN 
interaction EqnfT^lat / = 1/3 in a triangular lattice. (3) 
The transition from the SF to the CDW-VB phase in a 
Kagome lattice at / = 2/3 was studied by the DVM in 25 
and this paper in the easy plane limit v < 0, w > of the 
effective action Eqnllgl and the QMC on the EBHM of 
hard core bosons with NN interaction EqnJTTJat / = 2/3 
in a Kagome lattice. 

In this paper, by using the method developed in this 
paper, we confirmed the first two and established the 
third explicitly for the first time. The first connection is 
on Ising limit. It was known that the transition is likely 
to be first order in the Ising limit. The last two are the 
only two known examples so far realized in the easy plane 
limit where there is some chance to realize the DCQCP 22 
from the DVM. The last one is the only known example 
involving some VB order where there is a better chance 
to realize the DCQCP from the DVM. The advantage 
of the frustrated lattices over the bipartite lattices are 
that local VBS order can be realized in the former even 
without any ring exchange interaction. Unfortunately, 
both cases came up as negative. The first is a strong 
first order transition. From the DVM in Sect.II-B-1, it 
is related to a sharp change of the saddle point structure 
of the dual gauge field from the SF to the stripe CDW 
phase. The second is a very weak first order transition 
which needs some high power QMC technique such as 
double-peaked histogram to distinguish the very weak 
first order from a possible second order transition. From 
the DVM in Sect.III-B-1, we showed that the form of 
effective action in the CDW-VB side is identical to that 
in the SF side despite the reduced translation symmetry 
in the CDW-VB side. This gives an intuitive explanation 
why the transition is a very weak first order due to the 
very weak relevance of the w cos W term. 

We also inspired furthermore comparisons between the 
two approaches. For example, the CDW+VB order in 
Fig. 7 was not even searched in the QMC in 2 ^ in any 
parameter regimes. We suggest to measure the bond- 
bond correlation to identify this interesting state by the 
QMC in 23 . We also suggest the QMC simulations of 
the EBHM EqnfT2]with NNN interaction in a Kagome 




(a) (b) 



FIG. 18: The bond phase factors in bipartite lattices at half- 
filling, (a) dual square lattice (b) dual triangular lattice. 

lattice at 1/3 to search for a stripe solid phase and a 
stable stripe supersolid phase slightly away from 1/3. It 
is still not known how to realize the TVB phase in the 
triangular lattice Fig. 8 and in a Kagome lattice Fig. 10 
by a microscopic model. However, it was shown that the 
excitation spectra above the TVB phase is necessarily 
of a density wave origin, so there is no TVB supersolid 
phase in both triangular and Kagome lattice. We expect 
more examples along these comparisons are needed to 
have a complete picture of phases and phase transitions 
from the two complementary approaches. 
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Appendix A: Magnetic space group in a square 
lattice 

In this appendix, we apply the same method in the Sec- 
tion II to study the symmetry breaking insulating states 
in a square lattice. We only focus at the q — 2 case. 
We reproduced the results achieved ir>i using the oper- 
ator formalism. Obviously, our method leads to much 
more intuitive and physical identification of these states. 
Fig. 20a, b and Fig. 19a have also been previously identified 
in 3 ^ in the context of Z2 gauge theory of high temper- 
ature superconductors. However, we disagree with their 
arguments leading to Fig. 20b of the plaquette state. We 
believe this disagreement involves different physical inter- 
pretation of a local VBS state which are very important 
for the identifications of all the important phases in frus- 
trated lattices discussed in the main text. Furthermore, 
we also explicitly give elementary excitation spectra in 
the dimer and plaquette VBS state. 

The Harper's equation at q = 2 in a square lattice 3 - 
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FIG. 19: (a) The X-CDW in the Ising limit in a square 
lattice. The vortex fields ip(0, 0) = 1 - (y/2 - ip(0, 1) = 

1 + (V2- l)t,V(l,0) = (v^-l)-f,V(M) = (\/2-l)+i- 
The vortex current flowing counter clockwise is I — \/2 — 1 . 
All the bonds have the same strength K = v2 — 1. There 
are no frustrated bonds. The vortex current flowing counter 
clockwise is / = \/3/2. All the bonds have the same strength 
K — 1/2. There are no frustrated bonds. Following the 
current counting methods in the main text, one can see the 
density at the red site is n r = 1/2 + 4x1, that at the black 
spot is rib = 1/2 — 4x1. The x ~ n r — nt tunes the distance 
away from the SF to the CDW transition. When rib = 0, then 
n r — 1 which corresponds to the classical limit t/Vi = 0. (b) 
Ising limit in a honeycomb lattice. 



leads to c m (l = 0) = c m = (1, y/2- 1). uj = e l2lTf = -1. 
Denoting the dual square lattice site by x = (a±, 02). The 
vortex operator is: 

1 

m=0 
1 

= c m (-ir ai Ko+ei(-ir(-i) a2 ] (ai) 



m=0 



In the permutative representation: 

& = (0 o + 0i)/V2 
6 = -i{4>o-4>i)/V2 



FIG. 20: The ground states in the easy-plane limit v < in 
a square lattice (a) Dimer state at A > 0. The long arrow 
is \/2, the short arrow is y/2(y/2 — 1). The dashed line is 
the frustrated vortex bond with strength Kf v — —2(\/2 — 
l) 2 . The un-frustrated vertical vertex bond is K v = 2. The 
horizontal bond is Kh = 2(V2 — 1). The red bond which 
is perpendicular to the frustrated vortex bond is the boson 
valence bond in the direct lattice. The bosons are hopping 
back and forth along the red bond, (b) Plaquette state at A < 
0. The vortex fields ^(0,0) = 2[l+(v / 2-l)i], V(0, 1) = ^2[1+ 
(V2 - 0) = V2[l + (V2- l)i] = i>(0, 1), ^(1, 1) = 0. 

The arrows indicate both the magnitude and the phase of 
the vortex fields. Because the vortex field vanishes at (1,1), 
so there is a local superfluid ( or local VB order ) around 
this dual lattice point as indicated by the red plaquette. The 
kinetic energies emanating from the zero vortex field points 
are also zero. All the other bonds have strength K = 8(v / 2 — 
1). The bosons are hopping around the red plaquette. There 
are no frustrated bonds. There are no vortex currents in both 
(a) and (b) indicating no CDW order, so the boson densities 
are fixed at / = 1/2. 



(a) 



(b) 



(A2) 



FIG. 21: Local excitations in the dimer (a) and (b) plaquette 
states, the flipped bonds are denoted as green bonds. Com- 
pare with a local density excitation in the TVB in a triangular 
and Kagome lattice in Fig. 11. 



In the Ising case v > 0: <po = l,<f>x = 0, the vortex 
fields, kinetic energies and currents are shown in Fig. 19a 

In the easy plane limiti v < 0, there are also two cases 
depending on the signs of the quartic term A cos 49: 

(A) If A > 0, then O = ^ie m7r / 2 ,n = 0,1,2,3. The 
vortex fields, kinetic energies of the n = case with 4>q = 
01 = 1 is shown in Fig. 20a. It is a dimer state. 

(B) If A < 0, then O = 4> 1 e^ n+1 ^ 7T ^ ,n = 0,1,2,3. 
The vortex fields, kinetic energies of the n = case with 
4>i = 1 is shown in Fig. 20b which is plaquette state. 

The corresponding local low energy excitations in the 
dimer and plaquette states are shown in the Fig. 21. Be- 
cause the density remains uniform, so there is no associ- 
ated density wave excitation. The translational moving 



of these local excitations through the whole lattice leads 
to the excitations spectra shown in Fig. 24. 



Appendix B: Magnetic space group in a honeycomb 
lattice 

In this appendix, we apply the same method in the Sec- 
tion II to study the symmetry breaking insulating states 
in a honeycomb lattice whose dual lattice is a triangular 
lattice. We only focus at the q — 2 case. We reproduced 
the dimer state in the Fig. 23a achieved ir£. Unfortu- 
nately, the plaquette state in Fig. 23b was not identified 
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CO(q) 




CO(q) 




FIG. 22: The excitation spectra in (a) CDW. The lower 
branch of Ab appears due to the non-zero chemical potential 
fi in Fig. 3a. In the absence of this chemical potential, namely, 
(i — 0, then the Ab shown as a dashed line is well above the 
CDW mode A a . (b) CDW+ SS state. 



in£. We also work out the excitation spectra in all these 
insulating states and corresponding supersolid states. 

The Harper's equation at q = 2 in a honeycomb lattice 
leads to the vortex field in the Eqn.13 in 5 : 



■ik-x 



4>+{x) 



-ik-x 



«M£) 



(Bl) 



where k = (7r/3,7r/3). 

In the Ising case v > 0: 0+ = \,<j>- = 0, the vor- 
tex fields, kinetic energies and currents are shown in 
Fig. 19b. the corresponding excitation spectra are given 
in the Fig.22. 

In the easy plane limits v < 0, there are also two cases 
depending on the signs of the cubic term w cos 39: 

(A) If w > 0, then </>+ = (f>_e^ 2n+1 ^ 3 , n = 01)2 . 
The vortex fields, kinetic energies of the n = case with 
<p_ = 1 is shown in Fig. 23a. It is a dimer state. 

(B) If w < 0, then 0+ = 0_e l2mr / 3 ,n = 0,1,2. The 
vortex fields, kinetic energies of the n = case with (/>_ = 
1 is shown in Fig. 23b. It is plaquette state. 

Again, the corresponding local low energy excitations 
in the dimer and plaquette states in the honeycomb lat- 
tice can be similarly constructed as those in the square 
lattice Fig. 21. Because the density remains uniform, so 
there is no associated density wave excitation. The trans- 
lational moving of these local excitations through the 
whole lattice leads to the excitations spectra shown in 
Fig.22. From Eqn.22 in£, we can achieve the excitation 
spectra across the VBS to the VB-SS transition in Fig. 24. 
It is instructive to compare Fig. 24 with its counterpart 
across the CDW to CDW-SS transition Fig.22. 



Appendix C: Density operators on square latices 

In this appendix, we compare our gauge invariant den- 
sity operator, kinetic energy and current methods with 
the density operator formalism developed ini»2i. We 
pointed out some potential problems with the density op- 
erator formalism ini^ and also stress the different con- 
clusions reached by the two methods. 

In the dual vortex methods developed ini, how to char- 
acterize the CDW order of bosons on direct lattice in 





(a) 



(b) 



FIG. 23: Easy-plane limit in a honeycomb lattice (a) Dimer 
state. The dashed line is the frustrated vortex bond with 
strength Kf = — 1. The un-frustrated vertical vertex bond is 
K — 2. The blue bond which is perpendicular to the frus- 
trated vortex bond is the boson valence bond in the direct 
honeycomb lattice. The bosons are hopping along the blue 
bond, (b) Plaquette state. Because the vortex field vanishes 
at (1,0), the kinetic energies emanating from the zero vortex 
field points are also zero, so there is a local superfluid around 
this dual lattice point as indicated by the red honeycomb. 
The bosons are hopping around the red honeycomb. There 
are no frustrated bonds. All the other bonds have strength 
K = 3. There are no vortex currents in both (a) and (b) 
indicating no CDW order, so the boson densities are fixed at 
/ = 1/2. 
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q o 



(b) 



FIG. 24: The excitation spectra in (a) VBS. the VBS gap 
is given by the operator u>cos3# in Eqn.22 in 5 -. The lower 
branch of A appears due to the non-zero chemical potential 
/i in Fig.3b in 5 . In the absence of this chemical potential, 
namely, fi — 0, then the A is well above the VBS mode, (b) 
VB-SS state. 



terms of vortex operators on dual lattice is a little bit 
tricky. Ini, the density operator parameter was con- 
structed to be the most general gauge invariant bilinear 
combinations of the vortex fields: 



9-1 

r,s= — q 



i27T f (rx-\-sy) 



q-l 



Pre = S(\Q rs \cJ- rs ^Y,^^l+s^ lr (CI) 



where Q rs — 2irf(r,s) is the ordering wavevector and 
S'dQr-sl is a general form factor which can not be deter- 
mined from symmetry considerations and depend on the 
microscopic details. For simplicity without affecting the 
symmetry of the underlying state, it can be taken as a 
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FIG. 25: The possible bubble CDW phase in a triangular 
lattice in the easy plane limit v < 0,w > at / = 1/3. 
By counting the segments of currents along the bonds sur- 
rounding hexagons, paying special attentions to their counter- 
clockwise or clock wise directions, we can calculate the den- 
sities at the red and black points n r — 1/3 + 4x1 > 1/3, rib = 
1/3 — 2x1 < 1/3. The x ~ n r — rib > can be thought as a 
bubble CDW order parameter and can be tuned by the dis- 
tance away from the SF to the CDW transition in the Fig. 3a. 
When one tunes the density nj, = which stand for vacancies, 
then n r = 1 which stands for one boson. Note the vortex field 
at the center of the red loop is non- vanishing, only the current 
emanating from the center vanish due to the cancelations of 
currents from neighboring red sites, but the kinetic energies 
are positive. This indicates a local CDW order. 



a{)il){x),K v {x) = ip*(x + a 2 )e l27Tfai ip(x): 



( Pl K h ,K v )(x) 



Prs 



A 



K n = 



A = 



K v 



A"„ — 



9-1 



r,s—0 

9-1 



-lr 



1=0 



9-1 



1 * 



A h L ~ 

rs 



rn=0 

9-1 
1=0 

9-1 



-lr 



m(s+l) 



pr / j c m c m+r^ 
m=0 

9-1 

A; s [w-^- 1 >J2ti<t>i+ 

1=0 



-lr-\ 



1 



(C2) 



m=0 



Lorentzian S(Q) — qt+i ■ Then Eqn lCll was evaluated 
at the direct lattice point, the link points and the dual 
lattice points to represent the boson density, kinetic en- 
ergy and the ring exchange amplitudes respectively. As 
argued ir£, the link points of a square lattice is still a 
square lattice with a lattice constant y/2, the dual lattice 
points of a square lattice is still a square lattice withe 
same lattice constant. Putting the direct, link and dual 
lattices together forms a square lattice with a lattice con- 
stant 1/2. This may be the reason why the density oper- 
ator in Eqn lCll is the sum from — q to q — 1 and contains 
a factor lo~ ts I 2 which is non-periodic under / — > f + 1. 
As pointed out in£, it is very hard to extend the density 
operator defined ir>i in square lattice to other lattices. 
In this appendix, we give a simple way to construct the 
density and bond operators for any general q in terms 
of the dual vortex operator and then we will discuss the 
relations and difference between the two methods. 



The q eigenstates xi'J — 0, 1, • • • ,q — 1 which 
forms a q dimensional representation of the magnetic 
space group ( MSG ) [38] can be written as xi(^) = 
^J2 q m =o c rnU~ ml e i2 "^ mx+t ^ where x = (x,y). Then 

expand the vortex operator ip(x) = J21=o 4>i(x)xi{%) 
where 4>i{x) are the q order parameters. Just like 
Eqns l7l8l in the triangular lattices, one can construct 
the gauge invariant generalized density operator (GDO) 
p{x) — ip*(x)ip(x), the bond operators along the 
horizontal and vertical directions Kh{x) = ip*(x + 



By construction, all these operators transform like a den- 
sity, a bond or a current operator and is also periodic 
under /—>•/ + 1. We need evaluate these quantities only 
at dual lattice points to characterize the CDW and VBS 
orders in the direct lattice. 

Note that gauge invariant physical quantities con- 
structed in Eqns lC2l where the sum is from to q — 1 
and contain only / — > / + 1 periodic factors such as 
uj~ rs . However, in Eqn lCll the sum is from — q to q — 1 
and contains a factor uj~ rs / 2 which is non-periodic under 
/ -> /+ 1. Of course, the original EBHM in Eqn[T] is 
not periodic under / — > f + 1 , because all the interaction 
terms do depend on the filling factors /, however, we be- 
lieve the the symmetry of the insulating states should be 
periodic under / — > / + 1. This is what the EqnJCl] orig- 
inally designed for: characterize the symmetry breaking 
patterns of the insulating states only without the ability 
to spell out the microscopic details. As argued in£, al- 
though Eqn lCll seems work well in a square lattice and 
gives the same states shown in the appendix A, but it may 
not work in a honeycomb lattice. Similar density opera- 
tors in a triangular were constructed ir*2i, they gave the 
same CDW state in the Ising limit v > and the stripe 
CDW state in the easy plane limit v > 0, w < 0, but 
it leads to a bubble CDW state in the easy plane limit 
v > 0, w > in Fig. 2 5 which is completely different than 
the CDW-VB state shown in Fig. 7. This disagreement 
calls for some re-examinations of the density operators 
constructed iniSi. No such density operators were con- 
structed on honeycomb and Kagome lattices yet. 
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Appendix D: The vortex fields at a triangular 
lattice in the easy plane limit v < 0, w > 0. 

The vortex fields are shown graphically in the Fig. 5 
and have period 3 along both a\ and ai directions. Their 



specific values at x = (ai = 0, 1, 2, a-i = 0, 1, 2) are listed 
in the following two 3x3 matrices: 



2sin(7r/3)cos(7r/18)e" r / 9 
ip a (x)=\ 2sin(7r/3)cos(7r/18)e 427r / 9 
3sin(47r/9)e lTr / 6 

sin(7r/3)e- l7r / 6 
ifo(x) = I -sin(7r/3)e l7r / 6 



2 sin(vr/3) cos(7r/18)e-' i47r / 9 

2sin(7r/3)cos(7r/18)e~ l5,r / 9 
2 sin(7r/3) cos(7r/18)e^ l5,r / 9 -2 sin(7r/3) cos(7r/18)e l7r / 9 

-2i sin(7r/3) cos(2tt/9) 2 sin(7r/3) cos(2vr/9)e 47r / 6 

sin(7r/3)e 47r / 6 2 cos(27r/9)(sin(7r/9) + sin(27r/9))e l7r / 6 



-2icos(27r/9)(sin(47r/9) +sin(27r/9)) -2 sin(7r/3) cos(27r/9)e i7r / 6 2 cos(27r/9)(sin(47r/9) + sin(27r/9))e 47r / 6 



1 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. 
Sengupta, Phy. Rev. B 71, 144508 (2005). 

2 G. Murthy, D. Arovas, and A. Auerbach, Phys. Rev. B 55, 
3104 (1997). 

3 Longhua Jiang and Jinwu Ye, J. Phys. Cond. Matt. 18, 
6907 (2006). 

4 Jinwu Ye, cond-mat/0503113 , unpublished. 

5 Jinwu Ye, Nucl. Phys. B 805, 418 (2008). 

6 D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. 
Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

7 M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, T. 
Bloch, Nature 415, 39 (2002). 

8 G. Grynberg, B. Lounis, P. Verkerk, J.-Y. Courtois, and 
C. Salomon, Phys. Rev. Lett. 70, 2249 (1993). 

9 L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts , 
H. Fehrmannl, and M. Lewenstein, Phys. Rev. Lett. 93, 
030601 (2004); B. Damski, H. Fehrmann, H.-U. Everts, M. 
Baranov, L. Santos, and M. Lewenstein, Phys. Rev. A 72, 
053612 (2005). 

10 A. Griesmaier, J. Werner, S. Hensler, J. Stuhler and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005) 

11 K.-K. Ni, et ai, Science 322, 231 (2008). 

12 B. Capogrosso-Sansone, et ai, larXiv:0 906.2009 

13 E. Urban et ai, Nat. Phys, in press; A. Gaetan. et ai, Nat. 
Phys. in press. Jan. 2009 

14 G. Pupillo, et ai, larXiv:1001.0519l 

15 K. Baumann, C. Guerlin, F. Brennecke, T. Esslinger, 
larXiv:0912.326"Tl to appear in NATURE, 

16 H. P. Bchler, M. Hermele, S. D. Huber, M. P. A. Fisher, 
and P. Zoller, Phys. Rev. Lett. 95, 040402 (2005). 

17 R. G. Melko, et ai, Phys. Rev. Lett. 95, 127207 (2005). 

18 D. Heidarian and K. Damle, Phys. Rev. Lett. 95, 127206 
(2005). 

19 S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 
(2005). 

20 M. Boninsegni and N. Prokof'ev, Phys. Rev. Lett. 95, 
237204 (2005). 



21 A. A. Burkov and L. Balents, Phys. Rev. B 72, 134502 

(2005) . 

22 T. Senthil et.ai, Science 303, 1490 (2004). 

23 R. G. Melko, A. Del Maestro, A. A. Burkov, Phys. Rev. B 
74, 214517 (2006). 

24 Jing Yu Gan, Yu Chuan Wen, Jinwu Ye, Tao Li, Shi-Jie 
Yang, Yue Yu, Phys. Rev. B 75, 214509 (2007). 

25 K. Sengupta, S. V. Isakov, and Y. B. Kim, Phys. Rev. B 
73, 245103 (2006). 

26 S. V. Isakov, S. Wessel, R. G. Melko, K. Sengupta, and Y. 
B. Kim, Phys. Rev. Lett. 97, 147202 (2006). 

27 K. Damle and T. Senthil, Phys. Rev. Lett. 97, 067202 

(2006) . 

28 Jinwu Ye, larXiv:0804.3429l unpublished. 

29 Jinwu Ye, Phys. Rev. B 71, 125314 (2005). 

30 A. A. Burkov and E. Dernier, Phys. Rev. Lett. 96, 180406 

(2006) . 

31 C. Lannert, M.P.A. Fisher and T. Senthil, Phys. Rev. B 
63, 134510 (2001). 

32 Jinwu Ye and Longhua Jiang, Phys. Rev. Lett. 98, 236802 

(2007) ; Jinwu Ye, Phys. Rev. Lett. 97, 236803 (2006); 
Jinwu Ye, Annals of Physics, 323, 580 (2008). 

33 Jinwu Ye, J.M. Zhan g, W.M. Li u, Keye Zhang, Yan Li, 
Weiping Zhang, I arXiv: 1001. 32301 

34 F. Hbert et.ai, Phys. Rev. B 65, 014513 (2001). 

35 P. Sengupta, et.ai, Phys. Rev. Lett. 94, 207202 (2005). 

36 G. G. Batrouni, F. Hebert, R. T. Scalettar, Phys. Rev. 
Lett. 97, 087209 (2006); G.G. Batrouni, F.F. Assaad, R.T. 
Scalettar, P.J.H. Denteneer, Phys Rev A 72, 031601(R) 
(2005). 

For a review, see C. Lhuillier and G. Misguich, 
|arXiv:cond-mat/0109146| 
38 Jinwu Ye, Ph ys. Rev. Lett. 97, 125302 (2006); EPL 82, 
16001 (2008), |cond-mat/0603269| 



37 



Jinwu Ye, Jour, of Low Temp. Phys: 158, 882 (2010). 



